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Abstract 

Previous studies showed that for conditions simulating an 
aircraft encountering super-cooled water droplets the droplets 
may splash before freezing. Other surface effects dependent on 
the water surface tension may also influence the ice- accretion 
process. Consequently, the Weber number appears to be 
important in accurately scaling ice accretion. A scaling method 
which uses a constant- W eb er-numb er approach has been 
described previously; this study provides an evaluation of this 
scaling method. Tests are reported on cylinders of 2.5 to 15 -cm 
diameter and NACA 00 1 2 airfoils with chords of 1 8 to 53 cm in 
the NASA Lewis Icing Research Tunnel (IRT). The larger 
models were used to establish reference ice shapes, the scaling 
method was applied to determine appropriate scaled test 
conditions using the smaller models, and the ice shapes were 
compared. Icing conditions included warm glaze, hom glaze 
and mixed. The smallest size scaling attempted was 1/3, and 
scale and reference ice shapes for both cylinders and airfoils 
indicated that the constant- Weber-numb er scaling method was 
effective for the conditions tested. 

Nomenclature 

A c Accumulation parameter, dimensionless 

c Cylinder diameter or 2x airfoil leading-edge radius, cm 

c p Specific heat, cal/gm K 

h c Convective film heat-transfer coefficient, 

cal/sec m 2 K 

k Thermal conductivity, cal/sec m K 

K 0 Modified droplet inertia parameter, dimensionless 
LWC Liquid-water content, g/m 3 

M Mach number, dimensionless 

MVD Median volume droplet diameter, pm 

n Freezing fraction, dimensionless 

p Ambient static pressure, nt/m 2 

P w vapor pressure of water, nt/m 2 

R a Gas constant for air, 287.0 nt m/kg K 

Re Reynolds number, dimensionless 

t Temperature, °C 

T Absolute temperature, K 

V Airspeed, m/s 

We Weber number, dimensionless 

p 0 Droplet collection efficiency, dimensionless 

d Droplet median volume diameter, pm 

<p Droplet-energy transfer term in energy equation, K 
0 Air-energy transfer term in energy equation, K 

Ay Latent heat of vaporization, cal/gm 


p Density, gm/cm 3 

ji Viscosity, gm/cm s 

r Icing time, min 

Subscripts: 
a Air 

d Based on droplet size 

R Reference size and conditions 

st Static 

surf Surface 

S Scale size and conditions 

tot Total 

w Water 


Introduction 

Icing wind tunnels are a valuable tool to simulate natural icing 
encounters to obtain data quickly and inexpensively. In contrast 
to flight testing, icing tunnels permit the selection and control of 
temperature and cloud conditions. However, tests in icing 
tunnels may also require significant compromises. In particular, 
limitations on model size often require that scaled test hardware 
be used, and limitations on coverage of the full FAR25 
Appendix C may prevent testing at the desired liquid-water 
content, LWC , and median-volume droplet diameter, d. 

A number of scaling laws have been proposed over the past 40 
years in an effort to permit simulation of natural flight in icing 
tunnels. These laws, or methods, require that the scaled test 
conditions be chosen such that the values of several parameters 
important to ice accretion be matched to the respective full-size 
values. In this paper, the desired, or full-size conditions will be 
called the 'reference' conditions. Scaling laws have matched 
such parameters as the accumulation parameter,/!^ the modified 
droplet inertia parameter, K Qy and various terms included in the 
Messinger 1 heat-balance equation, such as the freezing fraction. 
Scaling laws developed by Charpin and Fasso 2 and by Ruff 3 
have been used extensively; however, they have been found to be 
less effective for glaze ice conditions than they are for rime 4,5 . 

In 1988, Bilanin 6 published a Buckingham-p-Theorem analysis 
of icing and showed that in addition to the parameters 
traditionally considered in the development of scaling laws, the 
scale Mach, Reynolds and Weber numbers should also be 
matched to their respective reference, or full-size, values. In 
1 995, Bilanin and Anderson 7 showed inferential experimental 
evidence that droplet splashing on the surface during glaze ice 
accretions may be important. In that paper, a scaling method 
was presented in which the Weber number, We y was matched in 



addition to the parameters used in the traditional methods. By 
maintaining the same We between reference and scaled 
conditions, surface-water effects and interactions between 
impinging droplets and surface water for the scale test match 
those in the reference situation. 

Reference 7 presented comparisons of ice shapes from tests 
performed with cylinders of 2.5- to 1 5-cm diameter using the 
constant- We scaling method. Scale ice shapes generally 
matched reference shapes for sizes as small as 1/3 the reference. 
However, in that study, experimental problems caused the test 
We for scale conditions to differ from the reference values by as 
much as 28%. In the study reported here, several of those tests 
with cylinders were repeated with the We properly matched. 

The droplet trajectory analysis used in the development of the 
scaling equations was derived by Langmuir and Blodgett 8 for 
cylinders. To verify that the scaling equations also apply to 
airfoils, additional tests were performed as part of the present 
study using NACA 0012 airfoils of 53.3-cm (21-in) , 26.7-cm 
(10.5 in) and 17.8-cm (7-in) chords. Tests with these models 
provided reference ice shapes (using the 5 3 -cm chord) and 1/2 
and 1/3 scaled results for comparison. Reference test conditions 
included air speeds of 67 and 89 m/s (150 and 200 mph), total 
temperatures of -1 1 to -4°C (1 2 - 24°F), water droplet sizes of 
30 and 40 mm median-volume diameter, and cloud liquid-water 
contents of .55 to 1 .0 g/m 3 . These conditions produced mixed 
and glaze ice. 

All tests were performed in the NASA Lewis Research Center 
Icing Research Tunnel (IRT). This paper compares the scale 
and reference ice shapes which resulted from these tests. 

Description of Experiment 

NASA Lewis Icing Research Tunnel. The IRT is shown in 
figure 1 . It has been described in reference 9. The IRT has a 
test section width of 2.7 m (9 ft) and a height of 1 .8 m (6 ft.) It 
is capable of operation at test-section airspeeds of 1 60 m/s (350 
mph) and higher, depending on model size. A refrigeration 
system permits accurate control of the test-section temperature 
from-40to4°C (-40 to 40° F.) A water-spray system 10 with 8 
spray bars provides the ability to control test-section 
liquid-water content from .2 to 3 g/m 3 and droplet median 
volume diameters from 1 5 to 40 mm. 

Two sets of spray nozzles, known as the mod- 1 and standard 
nozzles, are used in the IRT to provide different ranges of 
liquid- water content and droplet size 10 . For the present tests, 
only the mod-1 nozzles were used. 

Test Hardware. Ice accretion was measured on circular 
cylinders and NACA 00 1 2 airfoils mounted in the center of the 
IRT test section. Hollow aluminum cylinders with 15.2-, 7.6-, 

5. 1 - and 2. 5-cm (6-, 3-, 2- and 1 -in) outside diameters were 


used. Figure 2 shows how each cylinder was positioned in the 
IRT test section. A retractable shield was located in front of the 
cylinder to protect it from spray during the spray-bar start-up 
period when the water and air pressures were stabilizing. The 
shield was retracted rapidly into the tunnel ceiling by remote 
operation of a hydraulic actuator when the spray conditions had 
stabilized. 

The NACA 00 1 2 airfoils had chords of 53.3 cm (2 1 in), 26.7 cm * 
(10.5 in) and 17.8 cm (7 in). The 53.3-cm-chord airfoil was 
made of fiberglass while the 26.7- and 1 7.8-cm-chord airfoils 
were machined from solid aluminum. Figure 3 shows the 53.3- 
cm-chord, full-span model mounted on the IRT turntable. 

Because of the way this airfoil was installed, the retractable 
shield could not be used with it. 

The start-up shield was used with both of the smaller airfoils. 

Figure 4 shows the 26.7-cm-chord airfoil in the IRT test section 
with the shield. In Figure 4(a) the shield is lowered to protect 
the model from spray until cloud conditions became stable. 

Figure 4(b) shows the shield raised to accrete ice on the model. 

The installation for the 1 7.8-cm-chord airfoil was identical. 

Both of these aluminum models had spans of 1 22 cm (4 ft). 

All airfoils were mounted at 0° angle of attack for these tests. 

Test Procedure. Tests were performed by first establishing the 
desired velocity and temperature. Water spray was then 
initiated. When the shield was used, it was raised when the 
spray conditions had stabilized, and the spray timer was started 
at this time. When the shield was not used, .2 min were added to 
the required spray time to adjust for the period of spray transient. 
When the prescribed spray period was completed, the spray was 
shut off and the tunnel brought to idle to permit personnel entry 
into the test section to trace the ice shape. 

The ice shape was recorded by first melting a thin slice through 
the ice normal to the model span-wise axis. The shape was then 
traced onto a cardboard template; the shape was later digitized 
for computer storage to generate plots for comparisons. After the 
ice shape was recorded, the model was cleaned and the 
procedure repeated for the next spray condition. 

Spray Stability. For each of the airfoil tests, tunnel and cloud 
conditions were recorded at 2 -sec intervals. Recording started 
with the initiation of the icing cloud for tests in which the model 
shield was not used (i.e., tests with the 53.3-cm-chord airfoil 
only) and with the lifting of the shield in all other tests. These 
data provided a means to determine the time to stabilize the 
spray in the case of the non-shielded tests as well as a check of 
the stability of the icing spray for both shielded and non-shielded 
starts. 

Figure 5 shows the cloud liquid-water-content history for a 
typical test without the start-up shield. The LWC was calculated 
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from the average tunnel velocity and the air pressure and water- 
air differential pressure for each spray bar at each recording 
interval. The calculation used calibration equations based on 
liquid-water content measurements in the IRT test section. The 
result for each spray bar is shown with a different symbol. The 
average of the 8 LWCs increased rapidly to reach the set 
condition of .55 g/m 3 in about 30 sec. After an overshoot of 
about 3%, the LWC returned to the set condition within 2 min of 
* spray initiation. For tests with the shield, the typical range of 
LWC variation after the shield was raised was ± .4%. For all 
tests without the shield, desired spray times were increased by .2 
min to compensate for the low LWC during the start-up 
transient. The desired spray time for the test of figure 5 was 9.9 
min before adding the start-up transient allowance. When the 
average product of LWC and time from Figure 5 was integrated 
over the actual 10.1 -min spray time, it was found to be within 
2.5% of the value of the desired product of .55 g/m 3 and 9.9 min. 
The spray times reported here are those without the .2 min 
compensation. 

Figure 6 gives the history of the MVD for each of the 8 spray 
bars for the same 10-minute icing test without the start-up 
shield. As with the LWC calculation, equations used were from 
calibrations of the cloud using drop-size measurements in the 
IRT test section. About 2 min was needed to achieve the 
required 40pm drop size with good stability. The average drop 
size for the cloud climbed to 35pm in less than 30 sec, overshot 
to 45pm in less than 1 min and reached 40pm again after about 
2 min. Thus, for 95% of the spray time, the drop size was within 
± 5pm of the correct value. This result was typical of the 
unshielded tests reported here. For tests using the shield, the 
typical variation of MVD was less than ± . 1 pm once the shield 
was raised to expose the model to the icing cloud. 

Typical IRT velocity histories are shown in Figure 7. The 
velocity was computed from the total pressure and the total-static 
differential pressure measured with each of two pitot probes 
located on either side of the entrance to the IRT test section. 
Velocities for the unshielded test whose LWC and MVD were 
given above is presented in Figure 7(a). The different symbols 
show how one pitot probe gave approximately .7 m/s (1.5 mph) 
higher velocity readings than the other. The minimum and 
maximum velocities calculated for this test differed by less than 
2 m/s (5 mph). 

When the shield was used to protect the model from the icing 
cloud during spray stabilization, the test-section blockage was 
higher with the shield lowered than when it was raised. The 
change in blockage which occurred when the shield was raised 
required the tunnel operators to adjust the velocity. Figure 7(b) 
shows a typical velocity history for a shielded-model test. Time 
started with the raising of the shield. Again, the two pitot probes 
gave velocities differing by about .7 m/s (1.5 mph). For this test 
the desired velocity was 88 m/s (196 mph); this value was 
reached within about 30 sec of raising the shield. The velocity 
stability exhibited in Figure 7(a) was never achieved, and the 
speed oscillated about the correct value for the duration of the 


test. However, the range from the minimum indicated from the 
lower-reading probe to the maximum from the higher was less 
than 3 m/s (6 mph). Both the oscillation and the magnitude of 
fluctuation shown in Figure 7(b) was typical of velocity 
observations for other tests using the shield. 


Uncertainty in Parameter Values. All test conditions reported 
are averages. When the shield was used to protect the model 
from spray during start-up, the averages were obtained after the 
shield was raised. For tests without the shield the average values 
of velocity, droplet size and liquid-water content were obtained 
after the stabilization of cloud conditions. To establish the total 
uncertanty in the velocity, droplet size and liquid-water content, 
the pressure transducer accuracy was combined with observed 
fluctuations in the measurements during each run. The 
maximum uncertainty in velocity was ±2.7%, in MVD, ±1 3 .4%, 
and in LWC, ±2.4%. 

The tunnel temperature was the average of 1 1 thermocouples 
located on the turning vanes just upstream of the spray bars. 
These thermocouples have an uncertainty of about ± .6°C 
(±1.5°F). In addition, stratification in the tunnel introduced an 
additional uncertainty of about ± .8°C (±2°F). This can be seen 
in Figure 8 which gives typical histories of each of the 1 1 
thermocouples located at turning- vane 'D' which is just upstream 
of the IRT spray bars. In this study, the test section temperature 
was not independently determined, and the temperatures 
reported are the averages of the 1 1 comer-D’ thermocouples. 
Thus, although some mixing of the flow may reduce the 
stratification, the uncertainty in temperature may be as much as 
±1.4°C(±3.5°F). 


Scaling Equations 

Reference 6 proposed a scaling method in which the scale 
droplet trajectory, Weber number. We , accumulation parameter, 
A cy freezing fraction, n ? and droplet energy transfer term, <p , were 
matched with the respective reference values. This method is 
similar to the AEDC method derived by Ruff 3 with the added 
restriction that We must be matched. The consequence of this 
restriction is that only the scale model size, c s > can be chosen; 
the scale velocity, V s is determined by the We. The We accounts 
for surface effects, and it is assumed that both Re and M can be 
neglected. The complete derivation is given in reference 6, and 
the working equations will be repeated here. 

The following four equations are first solved simultaneously to 
give the scale airspeed, V s . , scale static temperature, T s , scale 
static pressure, p s , and scale drop size, S s 


V S = V R 
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The scale liquid-water content, LWC S> was found by equating the 
scale and reference freezing fraction. The equation that results 
after equating scale and reference collection efficiency, and 
scale and reference droplet energy transfer term, <j>, is 


6- K* Vjt 

LWC S = 

s *o*K* v s 


( 5 ) 


From the Messinger equation, Bis the air energy transfer term: 


The final scale parameter needed is the spray time, r. It was 
found by equating the scale and reference accumulation 
parameters: 


LWC r V r c s 
S ~ R LWC s V s c r 


( 8 ) 


Table I presents several reference conditions and the 
corresponding scaled conditions when equations (1) - (8) were 
used to apply constant- IFb scaling. Sizes were scaled to 1/2 and 
1/3 of the reference values. Reference conditions were chosen to 
represent a range of ice types from near rime to warm glaze. 

It was noted above that constant- We scaling permits only the 
choice of scale size, unlike traditional scaling methods which 
allow the choice of both scale size and airspeed. Table I 
illustrates some additional potential disadvantages of this 
method. First, the scale airspeed is always higher than the 
reference value. To simulate an icing encounter at 90 m/s (200 
mph) with a 1 /3-scale test requires a tunnel capable of operating 
at about 140 m/s (310 mph). Mach-number effects , which have 
been neglected in this scaling method, are probably important 
for speeds much higher than this; thus, tests at scales smaller 
than 1/3 are probably impractical unless very low airspeeds are 
to be simulated. Second, the high airspeeds produce high total 
temperatures which can approach and exceed the freezing point 
of water. When ice is accreted at near-freezing temperatures the 
adhesion of the ice to the surface is greatly reduced, and 
aerodynamic forces may affect the shape. These effects are not 
included in the constant- We analysis; therefore the method 
should not be used when temperatures approach freezing. 
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( 6 ) 


where r is the recovery factor, taken as .875 in this analysis, 

Pw svrf* s va P or pressure of water at the surface of the model 
(i.e., at r w ) and p w is the vapor pressure of water in the 
atmosphere (i.e., at the absolute static temperature, T). The 
vapor pressures used for this study were from Pruppacher and 
Klett 11 . 

The scale and reference convective heat transfer coefficients, h c s 
and h CfR , in equation (5) were evaluated using the following 
empirical expression from Gelder and Lewis 12 : 


k a 

1.244 — 
c 


/ \ .4774 

Vc p a 


( 7 ) 


Results 

The constant-Weber-number scaling method was tested in the 
IRT using the ice-accretion test techniques discussed earlier. 
Scaled conditions from Table I were used with cylinders and 
NACA 0012 airfoils. This section reports the ice shapes 
recorded from those tests. For each test the ice shape was 
recorded at 2 locations along the span of the model: at 91 cm 
(36 in) from the tunnel floor (i.e., at the tunnel centerline) and at 
1 12 cm (44 in) (i.e., at 20 cm (8 in) above the tunnel centerline.) 
The cloud LWC calibration 10 was based at the center of the test 
section and cloud uniformity depends on tunnel airspeed; 
therefore, even with valid scaling methods, conditions which 
scale at the tunnel centerline may not scale elsewhere if the cloud 
uniformity changes significantly from reference to scale 
conditions. Thus, shapes at both locations will be reported. The 
x-y coordinates of the scaled ice shapes were multiplied by the 
appropriate scale factor to permit direct comparison with the 
reference shapes. Average test conditions are given with each 
figure. 
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Cylinder Tests. In the tests reported in reference 7 a transducer 
calibration error produced cloud conditions which gave scale We 
different from the reference values. Three tests from that study 
were repeated here with properly-scaled conditions. Cylinder 
tests were also made using reference conditions not reported in 
reference 7. 

Figure 9 shows results for tests of case 1 in Table I. This was a 
* warm-glaze condition. The reference ice shape for a 1 5.2-cm- 
diameter cylinder is shown as a solid line, and the 1 /2-scale test 
on a 7.6-cm cylinder is shown as a dashed line. Results for a 
second test of the 1 /2-scale condition, indicated with a dotted 
line, are given to show repeatability of the ice shape. Figure 9(a) 
gives the ice shapes recorded at the tunnel centerline position, 
and Figure 9(b) at 20 cm (8 in) above the centerline. The scale 
ice shapes at both locations matched the reference shapes within 
the repeatability indicated. 

The conditions tested for Figure 1 0 resulted in ice shapes with 
horn spread varying along the cylinder span in a cyclic manner. 
The hom angle of the reference shape (solid line) can be seen to 
vary from the centerline to 20 cm above centerline. The hom 
angle for the 1 /2-scale test also varied and did not precisely 
coincide with the reference shape. The agreement at the 
centerline position was somewhat better than that at 20 cm (8 in) 
above the centerline. 

Figure 1 1 gives results for a third case repeated from the tests of 
reference 7. This is again a warm glaze condition similar to that 
shown in Figure 9; however, the reference drop size for this test 
was 58 pm and the reference cylinder diameter was half that of 
Figure 9. Figure 1 1(a) shows that the 2/3 scale result matched 
the reference quite well while the 1/3 scale shape gave only fair 
agreement at the tunnel center. No reference shape was 
recorded at 20 cm above the tunnel center, but it can be seen 
from Figure 11(b) that the 2/3 and 1/3 scale shapes appear to be 
similar. 

The conditions reported in Figure 12 differed from those of 
Figure 1 1 in that the reference droplet size was reduced from 58 
to 40 pm and the liquid-water content was reduced from .89 to .8 
g/m 3 . The reference model was, again, a 7.6-cm (3-in) cylinder, 
for this test, only 2/3-size scaling was possible because 1/3- 
scaling produced conditions not achievable in the IRT. The 
reference test was repeated (solid and dashed lines). The 
repeatability was fair at the centerline (figure 12(a)) and 
excellent at the higher position (figure 1 2(b)). The match of the 
2/3 -scaled result (dotted line) was also excellent at both 
locations. 

Finally, Figure 13 presents a 2/3-size scaling test with the 
reference temperature reduced by 2.8°C (5°F) from the tests of 
Figure 12. Again, a repeat test was made at the reference 
condition. As for the conditions of figure 1 2, repeatability at 20 
cm (8 in) above the tunnel center was better than at the center. 
The 2/3-scaled test gave a tunnel-centerline ice shape that 
agreed better with the repeated reference test results than with 


the original reference shape. At the 20-cm- (8-in-) above- 
centerline location the scaled ice shape was a close match of 
both reference shapes. 

Airfoil Tests. All airfoil tests used a 53.3-cm- (21 -in-) chord 
NACA 001 2 as the reference model. Scaled ice shapes were 
obtained from tests with 26.7-cm- (10.5-in-) and 17.8-cm- (7- 
in-) chord NACA 0012 models (1/2 and 1/3 scale). 

Figure 14 compares 1 /2-scale and 1 /3-scale ice shapes with the 
reference shapes for a mixed ice shape at 89 m/s (200 mph). 
Shapes from both of the scaled tests matched the reference 
shapes at the centerline (Fig. 14(a)) and at 20 cm above the 
centerline (Fig. 14(b)). 

When the conditions of Figure 14 were repeated but with a 
2.8°C (5°F) higher temperature, the horn-glaze condition shown 
in Figure 1 5 resulted. The hom-glaze ice shape is often difficult 
to reproduce in a scale test. The constant- We method produced 
fairly accurate scale shapes at the tunnel center (Fig. 15(a)). At 
20 cm (8 in) above the centerline, the 1 /2-scale test reproduced 
the reference fairly well, but the 1 /3-scale test appeared to 
accrete too much ice. 

An additional increase in temperature of 2.8°C (5°F) over the 
tests of Figure 15 gave the hom-glaze ice shape shown in Figure 
16. Once again, the 1 /2-scale test simulated the reference ice 
shape very accurately while the 1 /3-scale appeared to accreted 
too much ice at both locations on the model. 

The test results given in Figure 17 are for a hom-glaze condition 
at 67 m/s (150 mph). Both the 1/2- and 1/3 -scale tests gave 
shapes that agreed well with the reference at the tunnel center 
(Fig. 17(a)). At the 20-cm (8-in) position (Fig. 17(b)), both 
scale tests appeared to accrete somewhat too much ice. 

An increase in temperature of 2.8°C (5°F) over that shown in 
figure 17 gave the ice shapes of Figure 18. Figure 18(a) shows 
that 1/2 and 1/3 scaling both gave excellent agreement with the 
reference ice shape. At 20 cm (8 in) above the centerline (Fig 
1 8(b)), the 1 /2-scale test gave only fair agreement with the 
reference while the 1 /3-scale did not match. 

Although the scaled ice shapes did not match the reference in 
every test, there was sufficient agreement to provide strong 
evidence of the validity of the constant- We scaling method. 

Concluding Remarks 

Tests have been performed using both cylinders and NACA 
0012 airfoils of various sizes to evaluate the constant- We scaling 
method. All models were mounted vertically in the IRT test 
section. Reference sizes included 15.2-cm- (6-in-) and 7.6-cm- 
(3-in-) diameter cylinders and 53.3-cm- (21-in-) chord airfoils. 
Comparison tests with conditions determined by the constant-JFe 
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method were made with models scaled to 2/3, 1/2 and 1/3 the 
reference size. Reference conditions were chosen to produce a 
variety of accreted ice shapes including both rime and glaze ice. 
Ice-shape cross sections were recorded at two locations along 
the span of the model: at the center of the tunnel (9 1 .4 cm (36 
in) from the floor) and at 20 cm (8 in) above the center. For 
most tests, scaled ice shapes with 2/3- and 1/2 -size models 
closely matched the reference shapes. Somewhat less success 
was experienced with the 1/3 -size tests, although many of these 
also matched the reference ice shapes. 

As with many scaling techniques, the constant- Jfe method 
requires the matching of the scale and reference freezing 
fraction. It was shown in this study that the uncertainty in tunnel 
temperature, on which the freezing fraction depends, may be as 
much as ±1 .4°C (±4°F). It was shown in reference 5 that a 
change in temperature of this magnitude can have a significant 
effect on ice shape. Clearly, then, improvements in temperature 
measurement are needed for future scaling studies and for 
accurate testing using any scaling method. 

In general, the best match of ice shapes was found at the mid- 
span (tunnel center) position. The constant- ffe method requires 
increased airspeeds for reduced-size models, and the cloud 
uniformity in the IRT changes with airspeed 10 ; consequently, 
scale and reference ice shapes recorded at the 20-cm- (8-in-) 
above-midspan position matched less frequently than those from 
the mid-span position. Cloud uniformity could be a significant 
concern for scaling when it’s necessary to determine ice shapes 
on more than just the midspan of the scaled model. 

There is no definition at the present time of how closely ice 
shapes need to match to claim successful scaling. For this study, 
adequate scaling was assumed if scaled shapes agreed with 
reference shapes within the repeatability of those shapes. The 
evaluation was entirely subjective, however. In the absence of 
either a quantitative description of ice shapes or the 
identification of critical features of ice shapes, assessment of 
scaling methods will continue to be based on subjective 
judgement. 

The success of the constant- We scaling method suggests that 
surface effects, which may include droplet splashing on impact, 
play a significant part in the ice-accretion process. A better 
understanding of such effects is needed to incorporate into ice- 
accretion-prediction models as well as scaling methods. 

In reference 7 an excellent match of ice shapes was 
demonstrated for scaling tests even when the scale We was as 
much as 28% lower than the reference. When the scale We was 
carefully matched to the reference value for the present study, the 
agreement between scale and reference ice shapes was no better 
than shown in reference 7. Further study is needed to determine 
how much the constant- ffe conditions can be relaxed and still 
produce acceptable scaling. 
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Table I 

Constant-lFe Scaling Conditions 


chord chord 
or or 

model cyl cyl 
diam, diam, 
cm in 


1 reference |cylinder | 15.2 


1/2 scale cylinder 


2 reference 


1/2 scale 


3 | reference [cylinder 


3 scale cylinder 5,1 


1/3 scale cylinder 2,5 


reference 


2/3 scale | cylinder | 5.1 


2.5 


5 reference | cylinder | 7.6 


2/3 scale | cylinder 1 5.1 




airfoil 


1/2 scale I airfoil 


1/3 scale airfoil 


9 reference airfoil 


1/3 scale 

airfoil 

reference 

airfoil 


tton 


t tot* 

V, 

V, 

°C 

°F 

°F 

m/sec 

mph 


-4.4 20.0 24.0 67.1 150.0 


-7.11 -3.1 19.3 26.4 89.1 199.4 


200.0 


117.4 262.7 


-6.8 -3.8 


- 2.2 




-10.01 -5.4| 14.0| 22.21 95.5 213.61 


-6.81 -3.9 


6 

reference airfoil 

53.3 

21 


17.8 


53.3 21 


26.7 | 10.5 


17.8 


53.3 


26.7 10.5 


17.8 


21 


10.5 


53.3 21 


26.7 10.5 


-8.8 3.7 


-16.3 -6.8 2.6 


-12.2| -8.2 10.01 


-12.9 


-13.6 


-5.5 15.0 


-10.1 -3.2 


- 11.1 


-11.5 -7.6 11.3 


-11.8 -6.5 


-8.3 

-8.7 


-9.1| -3.8 15.7| 


27.2 


12.2 89. 


16.2 117.9 263.7 


19.8 138.6 310.11 


17.2 89.4 200.0 


263.6 


309.9 


22.2 89.4 200.0 


26.2 


263.5 

29.8 1 138.5| 

309.8 


m 


0.205 3.80 2.77 


2.94 


0.367 1.11 


0.275 1.68 4.92 


1.07 


196.7 


230.51 16.9 


150.0 


23.2 


25.2 103.0 230.5 


0.318 0.65 


mm 


































































































































































Figure 1. NASA Lewis Icing Research Tunnel (IRT). 



Figure 2. Test Cylinder and Shield Mounted in IRT. 


Figure 3. 53.3-cm-chord NACA 0012 Airfoil Installed in IRT 
Test Section. 


8 






LWC, g/m 



m 


ik* 


i 


■?m mm 








U 


(a) Start-up Shield Lowered to Protect Model During (b) Start-up Shield Raised to Expose Model to Spray. 

Icing Spray Stabilization. 

Figure 4. 26.7-cm-chord NACA 0012 Airfoil Installed in IRT Test Section. 




E ^ 

zi- 
p' 30 

S 25 

20 



Time, sec 

Time, sec 

Figure 5. Liquid-Water-Content History in IRT During 

Typical 10-min Spray Without the Start-up Shield. 

Figure 6. Droplet Size History in IRT During Typical 
10-min Spray Without the Start-up Shield. 









Temperature, 



Time, sec Time, sec 


(a) Model not Shielded During Spray Start. (b) Shielded Model. Shield Raised at 0 sec. 


Figure 7* Typical Velocity History. 



0 100 200 300 400 500 600 0 100 200 300 400 500 600 

Time, sec Time, sec 


(a) Total Temperature, -2.3°C (27.8°F). (b) Total Temperature, -10.6°C (12.9°F). 

Figure 8. Typical History of 1 1 Thermocouples at IRT Turning Vane ’D\ 
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(a) Ice Shapes at Centerline. (b) Ice Shapes at 8 in Above Centerline. 



cyl diam 


*tot 

** 

*lot 

V 

V 

MVD 

LWC 

r 


cm 

in 

°C 

°c 

°F 

°F 

m/s 

mph 

pm 

g/m 3 

min 

■ Reference 

15.2 

6 

-7 

-4 

20 

24 

67 

150 

40.6 

1.11 

42.8 

We Ml Scale 

7.6 

3 

-7 

-3 

19 

26 

89 

198 

23.1 

1.27 

14.4 

We 1/2 (repeat) 

7.6 

3 

-7 

-3 

19 

26 

89 

199 

23.1 

1.26 

14.4 


Figure 9. Constant-lFe Scaling Applied to Cylinders. Case 1 in Table I. 



(a) Ice Shapes at Centerline. (b) Ice Shapes at 8 in Above Centerline. 




cyl diam 


‘tot 

‘st 

‘to, 

V 

V 

MVD 

LWC 

T 

4 


cm in 

°c 

°c 

°F 

°F 

m/s 

mph 

pm 

g/m 3 

min 


■ Reference 

15.2 6 

-12 

-8 

10 

17 

90 

200 

40.7 

.88 

40.3 


We Ml Scale 

7.6 3 

-13 

-6 

9 

22 

119 

266 

23.6 

.97 

13.8 


Figure 10. Constant-lFe Scaling Applied to Cylinders. Case 2 in Table I. 
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(a) Ice Shapes at Centerline. (b) Ice Shapes at 8 in Above Centerline. 



cyl diam 

ht 

hot 

‘s, 

hot 

V 

V 

MVD 

LWC 

r 


cm in 

°c 

°c 

°F 

°F 

m/s 

mph 

pm 

g/m 3 

min 

— Reference 

7.6 3 

-7 

-4 

20 

24 

67 

151 

58.0 

.89 

26.7 

We 2/3 Scale 

5.1 2 

-7 

-4 

20 

25 

79 

177 

42.7 

.95 

14.2 

We 1/3 Scale 

2.5 1 

-7 

-2 

19 

28 

103 

230 

25.0 

1.03 

5.0 


Figure 11. Constant- IFe Scaling Applied to Cylinders. Case 3 in Table I. 



(a) Ice Shapes at Centerline. 


(b) Ice Shapes at 8 in Above Centerline. 



cyl diam 

hi 

hot 


hot 

V 

F 

MVD 

LWC 

Z 


cm 

in 

°c 

°C 

°F 

°F 

m/s 

mph 

pm 

g/m 3 

min 

Reference 

7.6 

3 

-9 

-7 

15 

19 

67 

150 

40.0 

.80 

26.8 

Reference (repeat) 

7.6 

3 

-9 

-7 

15 

19 

67 

150 

39.9 

.80 

26.8 

We 2/3 Scale 

5.1 

2 

-10 

-7 

14 

20 

79 

177 

28.7 

.88 

14.3 


Figure 12. Constant -We Scaling Applied to Cylinders. Case 4 in Table I. 
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(a) Ice Shapes at Centerline. 


(b) Ice Shapes at 8 in Above Centerline. 



cyl diam 

*st 

*tot 

(sf 

‘tot 

V 

V 

MVD 

LWC 

T 


cm in 

°C 

°C 

°F 

°F 

m/s 

mph 

pm 

g/m 3 

min 

Reference 

7.6 3 

-7 

-4 

20 

24 

67 

150 

40.5 

.80 

26.8 

Reference (repeat) 

7.6 3 

-7 

-4 

20 

24 

67 

150 

39.1 

.80 

26.8 

We 2/3 Scale 

5.1 2 

-7 

-4 

20 

25 

79 

177 

28.8 

.86 

14.4 


Figure 13. Constant - We Scaling Applied to Cylinders. Case 5 in Table I. 



(a) Ice Shapes at Centerline. 


(b) Ice Shapes at 8 in Above Centerline. 



chord 

‘s' 

‘tot 

*st 

‘to! 

F 

V 

MVD 

LWC 

r 


cm 

in 

°C 

°C 

°F 

°F 

m/s 

mph 

pm 

g/m 3 

min 

Reference 

53.3 

21 

-15 

-11 

6 

13 

90 

200 

40.4 

.55 

9.9 

We 1/2 Scale 

26.7 

10.5 

-15 

-8 

4 

17 

117 

263 

23.0 

.62 

3.4 

We 1/3 Scale 

17.8 

7 

-16 

-6 

4 

21 

138 

309 

16.6 

.64 

1.8 


Figure 14. Constant - We Scaling Applied to NACA 0012 Airfoils. Case 6 in Table I. 




(a) Ice Shapes at Centerline. (b) Ice Shapes at 8 in Above Centerline. 



chord 

‘s' 

hot 

‘s' 

hot 

V 

V 

MVD 

LWC 

T 


cm 

in 

°C 

°c 

°F 

°F 

m/s 

mph 

|im 

g/m 3 

min 

" Reference 

53.3 

21 

-12 

-8 

10 

17 

90 

200 

40.5 

.55 

9.9 

We Ml Scale 

26.7 

10.5 

-12 

-5 

10 

22 

118 

263 

23.0 

.61 

3.4 

We 1/3 Scale 

17.8 

7 

-13 

-4 

8 

26 

138 

310 

16.6 

.62 

1.9 


Figure 15. Constant-IFc Scaling Applied to NACA 0012 Airfoils. Case 7 in Table I. 



(a) Ice Shapes at Centerline. (b) Ice Shapes at 8 in Above Centerline. 



chord 


hot 

‘s' 

hot 

V 

V 

MVD 

LWC 

T 

Y 


cm 

in 

°C 

°C 

°F 

°F 

m/s 

mph 

pm 

g/m 3 

min 


Reference 

53.3 21 

-9 

-5 

16 

23 

89 

200 

40.6 

.55 

10.0 

y 

We M2 Scale 

26.7 

10.5 

-10 

-3 

14 

27 

118 

263 

23.0 

.59 

3.5 


We 1/3 Scale 

17.8 

7 

-10 

-1 

13 

30 

138 

310 

16.7 

.58 

2.0 



Figure 16. Constant-IFe Scaling Applied to NACA 0012 Airfoils. Case 8 in Table I. 
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(a) Ice Shapes at Centerline. (b) Ice Shapes at 8 in Above Centerline. 



chord 

ht 

* tot 

ht 

hot 

V 

V 

MVD 

LWC 

T 


cm 

in 

°c 

°c 

°F 

°F 

m/s 

mph 

pm 

g/m 3 

min 

— Reference 

53.3 

21 

-11 

-8 

13 

17 

67 

150 

41.0 

.65 

11.2 

We Ml Scale 

26.7 

10.5 

-11 

-8 

12 

19 

88 

196 

23.4 

.76 

3.7 

We 1/3 Scale 

17.8 

7 

-12 

-6 

11 

21 

103 

230 

16.9 

.81 

2.0 


Figure 17. Constant-^e Scaling Applied to NAC A 0012 Airfoils. Case 9 in Table I. 



(a) Ice Shapes at Centerline. (b) Ice Shapes at 8 in Above Centerline. 



chord 

ht 

^ tot 


hot 

K 

V 

MVD 

LWC 

T 


cm 

in 

°c 

°c 

°F 

°F 

m/s 

mph 

pm 

g/m 3 

min 

— Reference 

53.3 

21 

-8 

-6 

17 

21 

67 

150 

40.7 

.65 

11.2 

We 1/2 Scale 

26.7 

10.5 

-8 

-5 

17 

24 

88 

196 

23.2 

.74 

3.7 

We 1/3 Scale 

17.8 

7 

-9 

-4 

16 

26 

103 

231 

16.9 

.78 

2.0 


Figure 18. Constant-ire Scaling Applied to NACA 0012 Airfoils. Case 10 in Table I. 
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